UCT-TP-290/2012 
MZ-TH/12-32 
September 2012 



Corrections to the SU(3) x SU(3) 

(n: . . 

Gell-Mann-Oakes-Renner relation and chiral 

• n 

^ ! couplings Lg and i/J l 



(N 



oo 
O 
(N 



X 



J. Bordes*^"\ C. A. Dominguez P. Moodley J. Penarrocha*^''^ and K. 

Schilcher ('')'(^) 



% (")Departamento de Ffsica Teorica, Universitat de Valencia, and Institnto de Ffsica 

I . Corpuscular, Centro Mixto Universitat de Valencia-CSIC 

CD . 

r~| . ''Centre for Theoretical &: Mathematical Physics, University of Cape Town, Rondebosch 7700, 

South Africa 

^> \ ('^^Institut fiir Physik, Johannes Gutenberg-Universitat, Staudingerweg 7, D-55099 Mainz, 

■ Germany 



Abstract 

Next to leading order corrections to the SU{3) x SU{3) Gell-Mann-Oakes-Renner relation 
(GMOR) are obtained using weighted QCD Finite Energy Sum Rules (FESR) involving the 
pseudoscalar current correlator. Two types of integration kernels in the FESR are used to sup- 
^ ■ press the contribution of the kaon radial excitations to the hadronic spectral function, one with 

^ ■ local and the other with global constraints. The result for the pseudoscalar current correlator at 

zero momentum is tp5{0) = (2.8 it 0.3) x 10^'^ GeV^, leading to the chiral corrections to GMOR: 
= (55 lb 5)%. The resulting uncertainties are mostly due to variations in the upper limit 
of integration in the FESR, within the stability regions, and to a much lesser extent due to 
the uncertainties in the strong coupling and the strange quark mass. Higher order quark mass 
corrections, vacuum condensates, and the hadronic resonance sector play a negligible role in this 
determination. These results confirm an independent determination from chiral perturbation 
theory giving also very large corrections, i.e. roughly an order of magnitude larger than the 
corresponding corrections in chiral SU{2) x SU{2). Combining these results with our previous 
determination of the corrections to GMOR in chiral SU{2) x SU{2), 5-,^, we are able to deter- 
mine two low energy constants of chiral perturbation theory, i.e. Lg = (1.0 it 0.3) x 10^^, and 
H2 = —(4.7 ± 0.6) X 10~^, both at the scale of the /?-meson mass. 
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1 Introduction 



The light quark condensate, an order parameter of chiral-symmetry breaking, is known to play 
an important role in QCD as well as in chiral perturbation theory (CHPT) [Ij. In particular, 
the ratio of the strange to the up/down quark condensate, Rgq = {ss)/{qq), is a key QCD 
quantity measuring flavour SU{3) symmetry breaking in the vacuum [2]-[3], and entering in 
many QCD sum rule applications (for a review see ^j). This ratio is also related to two CHPT 
low energy constants, which in turn determine the next-to-leading order corrections to the Gell- 
Mann-Oakes-Renner (GMOR) relation [l],[6] 

^^5(0)11 = -{mi + m,){0\q-iqi+qjqj\0) = 2 fj, Uj, {1 - 5p) , (1) 

where i,j stand for (light) quark flavours, fp and Mp are respectively the decay constant and 
the mass of the pion or kaon, with fx/U = 1-197 ± 0.006, U = 92.21 ± 0.14 MeV [7j-[8j, and 
6p is the chiral correction. The relation between R^q and the CHPT low energy constants Lg 
and is M-W\ 

Rsq = = 1 + - 2^K - Mr, + ^{Ml - Ml){2Ll + HI) , (2) 

where < qq> \s the average of the up- and down quark condensates, Lg is a (physical) low-energy 
constant in the chiral Lagrangian to next-to-leading order, and 

IVlp IVlp 

^^ = 32^^"^' 

with the chiral renormalization scale. The constant Lg at a scale of the rho-meson mass has 
been estimated in chiral perturbation theory to next-to-leading order with the result [TT] 

Ll{v^ = Mp) = (0.88 ± 0.24) x 10"^ , (4) 

while a determination at order 0{p^) gives [I2j Ll{v^ = Mp) = (0.62 ± 0.20) x IQ-^ . The 
unphysical low energy constant H2 has been estimated in [TT] as 

Hl{v^ = Mp) = -(3.4 ± 1.5) X 10~^ . (5) 

These CHPT low energy constants are related to the next-to-leading order chiral corrections to 
GMOR as fohows 

Mt 

5^=A-l-{2Ll-Hl) and 5k = j^6^. (6) 

Jit -'"tt 

Using the results Eqs.(4)-(5) in Eq.(6) the following estimate was obtained in [llj 

6^ = (4.7 ± 1.7)% and 6k = (61 ± 22)% . (7) 

It is also possible to determine these corrections to GMOR entirely from QCD sum rules involving 
the correlator of the axial-vector current divergences at zero momentum. The feasibility of such 
an approach has been established long ago in a series of papers [I3]-[I8]. Recently, such an 



^In this respect there is a technical issue regarding logarithmic quark mass singularities which will be addressed 
in the next section. 
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independent determination of 5.,^ from weighted QCD Finite Energy Sum Rules (FESR), and in 
good agreement with the value above, has been obtained in |19j with the result 

= (6.0 ± 1.0)% . (8) 

Regarding 6k, the CHPT result in Eq.(6) already indicates that the chiral SU{3) x SU{3) 
corrections to GMOR are expected to be huge; in fact, much larger than typical chiral corrections 
such as the 10—15% deviation from the SU{3) xSU{3) Goldberger-Treiman relation [20]. Hence, 
it is important to confirm such a huge correction by determining it independently of the CHPT 
result in Eq.(6). In this paper we perform such a determination using modern weighted QCD 
FESR jlj, [H], [21] -[22] designed to reduce considerably the systematic uncertainties arising 
from a lack of direct experimental information on the pseudoscalar hadronic spectral function. 
In addition, by combining the result for 5k with our previous determination of 5-,^, Eq.(8), we 
are able to determine Lg and H2 independently, and entirely from QCD sum rules. 



2 Finite energy QCD sum rules 



The pseudoscalar current correlator V'5('?^)li) for quark flavours i = u and j = s, is defined as 

Mq'^Wu = i J d^xe^'^- < 0| r(j5(x)i5(0)) |0 > , (9) 

where < 0| is the physical vacuum and the pseudoscalar current density j5(x) is 

J5(a;) = (ms + mu) : s(x)i75u(x) : , (10) 

and rriu^s are the light quark masses. At zero momentum, and using a Ward identity one obtains 
the GMOR relation [T], [13]-^ 

^-5(0)1^ = -[ms{fi^) + m^{fi^)]{0\ss + uu\0){fi^) = 2 f], M], {1 -6k), (H) 

where /x is a renormalization scale. Using Wick's theorem in the time ordered product in Eq.(9) 
would lead to normal-ordered operators in the low energy theorem, Eq.(ll), and "05(0) would be 
a renormalization group invariant quantity. In principle, however, this is spoiled by the presence 
of mass-singular quartic terms in perturbative QCD(PQCD), as well as tadpole contributions, 
which would cancel out if the quark condensate becomes minimally subtracted instead of normal- 
ordered [ll]-[l8] (for a recent reassessment see 0]). In practice, though, one can safely neglect 
quartic mass corrections in the case of the up- and down-quark correlator, and consider normal- 
ordered quark condensates so that there is no residual //-dependence in '(/;5(0)|^. In the case of 
■05(0)1*, though, it would appear a priori that quartic strange-quark mass corrections could be 
sizable. However, a posteriori, and given the resulting uncertainties in the value of 6k it turns 
out that one can also disregard this issue which then becomes purely academic. 
The FESR are obtained by invoking Cauchy's theorem in the complex squared-energy plane. 
The pseudoscalar current correlator at zero momentum is then given by 

^5(0)A5(0) = i r M^ImMs)ds + ^ f Mfl^^^s)ds, (12) 
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where in the sequel we drop the quark-flavour indexes, Sth = M'^ is the hadronic threshold, the 
integration path C(|so|) is a circle of radius |so|, and A5(s) is an arbitrary analytic integration 
kernel to be described shortly. The radius |so| determines the onset of the hadronic continuum 
beyond the resonance region, which is assumed to be given by PQCD. Hence, the first integral 
in Eq.(12) involves the hadronic spectral function, and the second integral the QCD correlator. 
The latter is known up to five-loop order in PQCD |23j-|29j. Mass terms up to order m^, and 
leading order vacuum condensates with Wilson coefficients to two-loop order are also known [18J . 
A posteriori, these contributions turn out to be negligible on account of the overall uncertainty in 
ip5{0) (roughly at the 10% level), mostly due to variations inside the stability region in |so|; and 
to a lesser extent due to the uncertainty in the value of the strong coupling and the strange quark 
mass. The hadronic spectral function involves the kaon pole followed by its radial excitations, 
i.e. 

-Irmp^is) = 2flM^6is-Ml) + 1 Jm V5(s)|i?i?5 , (13) 
vr vr 

where the first two resonances, Ki = K{14:60), and K2 = i^(1830), have masses and widths 
known from experiment, i.e. Mj^^ = 1460 MeV, = 250 MeV, Mk^ = 1830 MeV, Tk^ = 
250 MeV [7j. This experimental information, though, is restricted to the values of the masses and 
the widths. Hence, a reconstruction of the resonance spectral function entirely from data is not 
feasible, leading to an unknown systematic uncertainty from model dependency. For instance, 
inelasticity and non-resonant background are realistically impossible to model. In an attempt to 
reduce this model dependency, it was first proposed long ago [30]- [33] to normalize pseudoscalar 
spectral functions at threshold using CHPT as a constraint. In the case of the kaonic channel, 
the CHPT constrained resonance spectral function, which includes the important resonant sub- 
channel A'* (892) — vr, is given by [l33J (there is a misprint in [33] which is corrected below) 

i Im Us)\.„ = § ^ jp|jL_ eis -Mi). (14) 



where 



/(.)= / ^(u-Ml)is-u)liMj, 



Ml 



u 



and 



= {Ml,-uY + Ml.Tl. ■ ^''^ 

The pion mass has been neglected above, in line with the approximation = to be made 
in the QCD sector, and the resonant sub-channel yo(770) — K is numerically negligible. The 
complete hadronic spectral function is then 

- Im Ms)\hao = 2flMi 5{s -MI)+^- Im Us)\kJ^^^^^^^^^^^ 
vr vr (i + Aj 

+ -lTa.i)^{s)\pQCD 0{s - So) , (17) 
vr 
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where Im tp5{s)\QCD is the perturbative QCD spectral function modehng the continuum which 
starts at some threshold sq, -BW^i,2(s) are Breit-Wigner forms for the two kaon radial excitations, 
normalized to unity at threshold, and A controls the relative importance of the second radial 
excitation. The choice A ~ 1 results in a reasonable (smaller) weight of the K(1830) relative to 
the K(1460). While further embellishments to this parameterization could still be made, they 
would not hide its strong model dependency and its related unknown systematic uncertainty. 
This problem has recently been solved by a judicious choice of the integration kernel, A5(s), in 
the FESR [1], [19], [21]-[22]. In fact, choosing a second degree polynomial constrained to vanish 
at the peaks of the two radial excitations makes this contribution to the FESR completely 
negligible in comparison with the uncertainties due to the value of Og, the strange quark mass, 
and variation of the results in the stability region in | sq I • In the kaon channel the expression for 
A5(s) is 

A5(s) = 1-aos-ai (18) 

with A^iMf) = A5(M|) = 0, which fixes ao = 0.768 GeV"^ and ai = -0.140 GeV""^. Other 
choices of kernels with local constraints do not lead to any appreciable improvement in the 
quenching of the resonant contribution to the FESR. The optimal choice of integration kernel 
is most likely application dependent. For instance, pinched kernels in FESR for the vector and 
axial- vector correlators [34] -[36] improve considerably the agreement between data and the first 
two Weinberg sum rules. This circumstance, often referred to as duality violation, is a con- 
tentious issue as some model analyses claim its presence [37]- [38], and other model independent 
analyses claim the opposite [39]. In the present case, a combination of Eq. pH]) and pinched 
kernels does not work as well as Eq. ()18p on its own ^22j. However, it is possible to achieve some 
degree of pinching using other type of kernels as described next. 



A different kernel with a global, rather than local constraint, has been used recently in connection 
with the corrections to GMOR in chiral SU(2), and it involves Legendre type polynomials [19j 
as explained next. We consider the general polynomial 



Pn{s) = 



(19) 



m=0 



with coefficients fixed by (i) imposing the normalization condition at threshold 

Pn (s = M|.) = 1 , (20) 
and (ii) by requiring that the polynomial Pn{s) minimizes the contribution of the continuum in 



the range 



(fe+i) 
4h ,so 



in a least square sense, i.e., 

s"" Pnis) ds = form = 0, , 



so 



(fc+i) 

th 



. n 



(21) 



The meaning of the threshold is as follows. The hadronic resonance spectral function can 

be split into the sum of a term involving the known kaon radial excitations, and a term involving 
unknown resonances, non-resonant background, etc., i.e. 



/mVf^(s) = ^/mVf (s) + ^(s-sJJ+^^) /mV'f^^^ls) 



(22) 



i=0 
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where the upper index k stands for the number of resonances exphcitly included in the sum. The 
polynomials obtained in this way arc closely related to the Lcgcndrc polynomials, so that Eq.(21) 
is exact on account of orthogonality. Specifically, denoting by Vn{p) the Legendre polynomials, 
their relation to the Pn{s) above is 



where the variable x{s) is 



- Kwh ' '''' 

defined in the range x G [— Ijl] for s G [■^th"''^^' ^o]' '^^^ normalization of the Vn{x) is such 
that Vi{x) = X, V2{x) = {f)x"' - 3a;)/2 , etc. The use of these polynomials leads to results for 
■i/'5(0) in a region of large values of sq {sq ~ 20 — 50 GeV^), where it is stable, and where PQCD 
is unquestionably valid. Since the Legendre polynomials exhibit zeros, this type of integration 
kernel provides some amount of pinching. For instance, the sixth-degree polynomial vanishes at 
^Jsq c± 2.0 GeV, 9.0 GeV, and other four places at higher energies. 
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Figure 1: Results for ^^5(0) in units of 10 ^ GeV^ as a function of so, for as(M^) = 0.326 
(A = 350 MeV), and using the local integration kernel, Eq.(18). 
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Figure 2: Results for ^5(0) in units of 10 ^ GeV^ as a function of sq, for as{M^) = 0.326 
(A = 350 MeV). Curves (a), (b), (c) correspond, respectively, to Legendre-type integration 
kernels P4, P5, and Pq in Eq.(19). 



3 Results and conclusions 



We consider first the FESR involving the locally constrained integration kernel, Eq.(18), and 
in Fixed Order Perturbation Theory (FOPT). In this case the strong running coupling and the 
quark masses at a scale |s| = |so| are considered fixed, thus only logarithmic terms contribute to 
the integrals. The renormalization group summation of leading logarithms is carried out after 
integration by setting the scale to /i'^ = — |so|, and all integrals can be done analytically. For the 
strong coupling we use the wide range c5J(Mt-) = 0.31 — 0.35, which covers many different results 
as reported recently in [40] (including the PDG average [7J). Using this range of values of Og 
to determine the strange-quark mass in a similar framework as here ^21] gives nT^\j^{2 GeV) = 
100 — 104 MeV, in the MS scheme (for recent reviews see [ll]-[42]). Alternatively, one could 
use results from lattice QCD (for a compilation see [7]) most of which are in good agreement 
within errors with [21]. Combining both methods we consider the wide range mj|jg^(2GeV) = 
94 — 104 MeV. In Fig. 1 we show a typical result for 1^5(0) as a function of |so| in the wide 
stability region |so| = 3-6 GeV^ corresponding to al(M^) = 0.326 {Aqcd = 350 MeV), and 
rn^l j^(2 GeV) = 100 MeV. Taking into account the uncertainties in and leads to 

^-5(0) = (2.8 ± 0.3) X 10-^ GeV^ , (25) 

which translates into the chiral corrections to GMOR 

6k = (55 ± 5)% . (26) 
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Proceeding to the second choice of integration kernel, i.e. the Legendre type polynomials, 
typical results for "05(0) ^-^e shown in Fig. 2 for the same input values of as and mj. Sta- 
bility is achieved in a very wide region already for a sixth-order polynomial with the result: 
■05(0) = (2.50 — 3.00) X 10~^ GeV^ , in excellent agreement with the FOPT determination, 
Eq.(25). The chiral correction, Eq.(26), confirms the CHPT determination of 

We discuss in closing a determination of the two CHPT constants Lg, and H2 at the scale of 
the p-meson mass. It was first pointed out in [43], and later reiterated in [2] and [4], that 
in order to obtain the ratio Rgq = {ss)/{qq) with some reasonable accuracy one should avoid 
deriving it from the ratio of V'5(0)m/''/'5(0)^, as this leads to a huge uncertainty. Instead, one 
should first determine the ratio of the scalar to the pseudoscalar correlator at zero momentum, 
RvA = V'(0)n/V'5(0)u5 and then use Rgq ^ {1 + RvA)/i^ — Rva)- The most recent determination 
in a framework similar to the one used here, i.e. using weighted QCD FESR, gives [1] 

Rsq = T^ = 0.6 ± 0.1 . (27) 

Using this value in Eq.(2) determines 2U^ + H2, and from the first relation in Eq.(6) one obtains 
2Lg — H2. Using the result for S-^ in Eq.(8), obtained in the same framework as here, gives 

LI = (1.0 ± 0.3) X 10"^ (28) 

HI = -(4.7 ± 0.6) X 10"^ , (29) 

in very good agreement with [11]. The value of Lg above is also in agreement within errors 
with recent determinations from lattice QCD (for a review see [S]): -^>g = (0.62 it 0.04) x 10~^ 
[35], = (0.85 ± 0.04) x lO'^ [36], and = (0.88 ± 0.08) x 10"^ [17], as well as with recent 
analytical determinations: Lg = (0.6 ± 0.4) x lO'^ [IS], and Lg = (0.64 ± 0.16) x lO'^ 
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